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Abstract 

We model the two-point correlation function of galaxies in a cold dark matter Universe 
by combining two powerful theoretical tools — dissipationless N-body simulations of dark 
matter clustering (specifically the GIF simulations carried out by MPIA and the Virgo 
Consortium) and semi-analytic modelling of galaxy formation. We construct catalogues 
of galaxies containing a wide range of information for each galaxy, including magnitudes 
Q^ ■ in various bands, star formation rates, disk and bulge sizes (from the semi-analytic model 

described by Cole et al [||) and, importantly, spatial positions and peculiar velocities from 
the N-body simulation. We then use this information to study the clustering properties 
of galaxies. 
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Whilst semi-analytic models have proven useful in modelling and predicting numerous proper- 
ties of the galaxy population (Kauffmann, White & Guiderdoni ||; Cole et al ||; Somerville 
& Primack [|TjJ they provide only limited information on the spatial distribution of galaxies, 
and so it has been difficult to study in detail the clustering properties of galaxies using these 
models. Recently, semi-analytic models have been combined with N-body simulations to pro- 
vide the required spatial information (Kauffmann, Nusser & Steinmetz ||10|| ; Governato et al ||; 
Kauffmann et al [|ll|]). We have used a similar technique to study the clustering of galaxies in 
CDM universes within well constrained semi-analytic models. The main difference between our 
technique and that of Kauffmann et al is that whilst they extract the merging history of 
each dark matter halo directly from the simulation, we construct this history using the extended 
Press-Schechter theory. We find that this gives the same statistical results as the Kauffmann 
et al method and allows us to resolve merger trees to much smaller masses. 

To construct a catalogue of galaxies containing spatial information we use the following 
procedure: (i) take the output from a dissipationless N-body simulation of dark matter and 
use a group finding algorithm (here we use the friends-of-friends algorithm with the standard 
linking length of b = 0.2) to locate bound, virialised haloes of dark matter of 10 or more 
particles (such groups have been shown to be stable by Kauffmann et al |ll|); (ii) determine 
the mass of each group, the position and velocity of its centre of mass and the positions and 
velocities of randomly selected particles within the group; (iii) for each group use a semi-analytic 
model of galaxy formation constrained to match the local B and K-band luminosity functions 
to determine the population of galaxies living within the dark matter halo; (iv) attach the 
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Figure 1: The local B and K-band luminosity functions from our model compared to vari- 
ous observational determinations. Note the good agreement between model (solid line) and 
observations (symbols) at the bright end. Our galaxy catalogues become incomplete at the 
faint end of the luminosity functions (faintwards of Mb — 5 log ft. ~ —17.5 in the B-band and 
Mk — 5 log ft ~ —20.5 in the K-band) due to the limited resolution of the N-body simulation. 
Our catalogues are constructed using only galaxies for which a complete sample is available 
within the model. 

central galaxy of the halo to the centre of mass of the group and attach any satellite galaxies 
to randomly selected particles within the halo so that galaxies trace mass within a given dark 
matter halo (which may not be exactly true in reality because of processes such as dynamical 
friction). 

This results in a galaxy catalogue which can be analysed to determine the clustering prop- 
erties of galaxies of any given luminosity, morphology, colour and so on. In fact by using this 
technique it is possible to produce catalogues of galaxies complete with spatial information (or 
alternatively redshifts, and angular coordinates) with any observationally motivated selection 
criteria. Furthermore, by identifying dark matter halos on the past lightcone of an observer a 
full kock galaxy redshift survey can be constructed. The luminosity functions determined from 
the galaxy catalogue are shown in Figure [l|. Although they become incomplete at the faint end 
because of the limited resolution of the N-body simulation there is good agreement between 
the bright ends and the observed luminosity functions. Since we only consider the clustering 
of galaxies for which our catalogue is complete the resolution limit is not important. We find 
however that an accurate match to the bright end of the luminosity function is important as 
it strongly constrains the resulting two-point correlation function. An example of the informa- 
tion produced by our model is given in Figure 0, which shows a slice through a ACDM N-body 
simulation upon which the positions of galaxies brighter than Mb — 5 log ft = —19.5 (we define 
the Hubble constant to be H = 100ft km s _1 Mpc -1 ) have been overlaid as circles. It can be 
seen quite clearly that the galaxies trace the mass to some extent. To determine exactly how 
well galaxies trace the underlying dark matter we estimate the two-point correlation function 
of these galaxies. 
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Figure 2: A slice through a ACDM N-body simulation. The volume shown is 141 x 141 x 
8 h~ x Mpc. Dark matter is shown by the greyscale, with the darker areas being the most 
dense. Overlaid are the positions of all galaxies brighter than M B — 51og/t = —19.5 indicated 
by circles. 



2 The two-point correlation function 

Shown in Figure [| is the two-point correlation function of galaxies brighter than Mb — 5 log h = 
-19.5 in our ACDM model (which has fi = 0.3, A = 0.7, h = 0.7 and a$ = 0.9). This is 
compared to the correlation function of the underlying dark matter and to the observationally 
determined correlation function of galaxies in the APM survey from Baugh Q. Note that all 
of the curves shown here are real space correlation functions. The galaxy correlation function 
shows many interesting properties. Firstly it is biased with respect to the mass correlation 
function, and furthermore this bias is scale dependent. On small scales there is in fact an 
antibias, which is exactly what is needed to reconcile the theory with the observed galaxy 
correlation function. The observed and model galaxy correlation functions agree over a wide 
range of scales, both showing approximate power law behaviour. The scale-dependant bias seen 
in these models arises from a complex interplay of effects. On large scales the bias is due to the 
intrinsic bias of dark matter halos in a CDM universe as described by Mo & White ||13|| . On 



smaller scales the bias is controlled by the way halos are populated with galaxies, specifically 
the variations in number of galaxies per halo for halos of a given mass. The Lagrangian radius 
exclusion of halos also affects the small scale bias. These issues are explored in greater detail 
by Benson et al 0, who also study f2o = 1 models and find that such models fail to reproduce 
the observed clustering of galaxies. 
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Figure 3: The correlation function of galaxies brighter than Mb — 5 log/i = —19.5 in a ACDM 
cosmology. Points with errorbars show the APM galaxy correlation function of Baugh Q . The 
dotted line shows the dark matter correlation function whilst the solid line shows the correlation 
function of galaxies in our model (the dashed lines to either side indicate the Poisson errors on 
this correlation function). 
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